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Abstract: This paper considers a cooperative OFDMA-based cog- 
nitive radio network where the primary system leases some of its 
subchannels to the secondary system for a fraction of time in ex- 
change for the secondary users (SUs) assisting the transmission of 
primary users (PUs) as relays. Our aim is to determine the cooper- 
ation strategies among the primary and secondary systems so as to 
maximize the sum-rate of SUs while maintaining quality-of-service 
(QoS) requirements of PUs. We formulate a joint optimization 
problem of PU transmission mode selection, SU (or relay) selection, 
subcarrier assignment, power control, and time allocation. By ap- 
plying dual method, this mixed integer programming problem is 
decomposed into parallel per-subcarrier subproblems, with each 
determining the cooperation strategy between one PU and one SU. 
We show that, on each leased subcarrier, the optimal strategy is to 
let a SU exclusively act as a relay or transmit for itself. This result 
is fundamentally different from the conventional spectrum leasing 
in single-channel systems where a SU must transmit a fraction of 
time for itself if it helps the PU's transmission. We then propose a 
subgradient-based algorithm to find the asymptotically optimal so- 
lution to the primal problem in polynomial time. Simulation results 
demonstrate that the proposed algorithm can significantly enhance 
the network performance. 

Index Terms: Cooperative communications, cognitive radio net- 
works, orthogonal frequency-division multiple-access (OFDMA), 
resource allocation, two-way relaying. 



I. Introduction 

Cognitive radio (CR), with its ability to sense unused fre- 
quency bands and adaptively adjust transmission parameters, 
has recently attracted considerable interest for solving the spec- 
trum scarcity problem |[T1|2|. A key concept in cognitive radio 
networks (CRNs) is opportunistic or dynamic spectrum access, 
which allows secondary users (SUs) to opportunistically access 
the bands licensed to primary users (PUs). Most of the works 
on dynamic spectrum access regard the secondary transmission 
as harmful interference and hence the SUs do not participate in 
the primary transmission. Recently, a new cooperation strategy 
between the primary system and the secondary system was pro- 
posed in and further investigated in [4|. Therein, the PU link 
leases its channel to the SUs for a fraction of time to transmit 
secondary traffic in exchange for the SUs acting as relays to as- 
sist the transmission of primary traffic. The spectrum-leasing 
based cooperation can improve the performance of both the pri- 
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mary and secondary systems and result in a "win- win" situation. 

The early works O and (4) on spectrum leasing only inves- 
tigated the time slot allocation in the single-PU, multi-SU, and 
single-channel scenario. More specifically, in [3|, one PU tar- 
gets at maximizing its rate while multiple SUs compete with 
each other to access the single channel. However, this scheme 
may result in an extreme case that the PU is aggressive and the 
SUs have no opportunity to access the channel. Recall that in 
CRNs, PUs are willing to share the spectrum resource with SUs 
if their quality-of-service (QoS) requirements are satisfied (HE). 
In (4), the PU maximizes its utility in terms of rate and revenue 
while the SUs competitively make decisions based on their rates 
and payments. Nevertheless, the virtual payment and revenue 
may lead to another extreme case that the PU provides all of the 
transmission time to the SUs on the single channel, which is not 
practical in CRNs. 

In this paper, we consider the general spectrum leasing 
and resource allocation problem in multi-channel multi-user 
CRN based on orthogonal frequency-division multiple-access 
(OFDMA). The motivation of using OFDMA is two-fold. First, 
OFDMA is not only adopted in many current and next gener- 
ation wireless standards but also a strong candidate for CRNs 
. The second is that OFDMA-based systems can flexibly in- 
corporate dynamic resource allocations in CRNs (e.g., (SHU). 
The primary system consists of multiple user pairs conducting 
bidirectional communication. The secondary system is a cellu- 
lar network consisting of a base station (BS) and a set of SUs. 
The two systems operate in a cognitive and cooperative man- 
ner by allowing the SUs to occupy certain subcarriers given that 
the QoS of the PUs are satisfied with the assistance of SUs as 
cooperative relays. 

As CRNs are typically hierarchical and heterogeneous, it is 
intuitive that if SUs can aggressively help PUs' transmission, 
then less subcarriers will be needed by the PUs to satisfy their 
QoS requirements, and as a result the SUs can access more sub- 
carriers for maximizing their own data rates. Meanwhile, as the 
communication in the primary system is bidirectional, the coop- 
eration of SU as relays can also bring network coding gain in 
the form of two-way relaying. Thus, more subcarriers can be 
leased to the SUs. The increased spectral efficiency is in turn 
transformed into cooperation opportunities. Optimizing such 
cooperative CRN has unique attractiveness and challenges as 
follows. 

Firstly, for the primary system, when relaying is necessary, it 
has to decide which cooperative transmission modes (one-way 
relaying and two-way relaying) to select and which set of SUs to 

1 ln two-way relay systems, a pair of nodes exchange information with the 
help of a relay node using physical layer network coding f9lfTTI. Two-way 
relaying can achieve much higher spectral efficiency than the traditional one- 
way relaying. 



choose, since it has higher priority in a CRN. Secondly, for the 
secondary system, it needs to schedule appropriate SUs to uti- 
lize the leased subcarriers for maximizing its total throughput. 
Moreover, for those SUs that not only be selected as relays but 
also be scheduled to transmit for themselves, the secondary sys- 
tem needs to balance their resource utilization. Thirdly, from the 
common perspective of the primary and secondary systems, it is 
crucial to determine which set of subcarriers to cooperate on to- 
gether with how much power and time slots to transmit signals, 
in order to satisfy the QoS requirements of the primary system. 

The main contributions of this paper are summarized as fol- 
lows: 

1. We propose an optimization framework for joint bidirec- 
tional transmission mode selection, SU selection, subcarrier as- 
signment, power control, and time slot allocation in the coop- 
erative CRNs. The objective is to maximize the sum-rate of all 
SUs while satisfying the individual rate requirement for each of 
the PUs. There are three distinct features in our optimization 
framework. First, by subcarrier assignment and allocating time 
slot between PUs and SUs in cooperation sessions, multiuser 
diversity can be achieved in both frequency domain and time 
domain. Second, as the communication in the primary system is 
bidirectional, we can exploit network coding gain in the form of 
two-way relaying to improve spectral efficiency via the SUs' as- 
sistance. Third, using the OFDMA-based relaying architecture, 
each PU pair can conduct the bidirectional communication by 
multiple transmission modes, namely direct transmission, one- 
and two-way relaying, each of them can take place on a different 
set of subcarriers. 

2. We show that in the multi-channel cooperative CRNs, the op- 
timal strategy is to let a SU exclusively act as a relay for a PU 
or transmit data for itself on a cooperated channel. This result 
fundamentally differs from the conventional cooperation in the 
single-channel scenario where a SU must transmit a fraction of 
time for itself if it forwards the PU's transmission on the chan- 
nel. 

3 . Using the Lagrange dual decomposition method, the joint op- 
timization problem is decomposed into parallel per-subcarrier- 
based subproblems. An efficient algorithm is proposed to find 
the asymptotically optimal solution in polynomial time. 

The remainder of this paper is organized as follows. Sec- 
tion II introduces the optimization framework, including system 
model and problem formulation. Section III presents the de- 
tails of the Lagrange dual decomposition method for the joint 
resource-allocation problem. Section IV provides the simula- 
tion results. Finally, we conclude this paper in Section V. 

II. Optimization Framework 

We consider an OFDMA-based CRN where the primary sys- 
tem coexists with the secondary system as shown in Fig.[T] The 
primary system is an ad hoc network, consisting of multiple user 
pairs with each user pair conducting bidirectional communica- 
tions. The secondary system of interest is the uplink of a single- 
cell network where a BS communicates a set of SUs. Note 
that the downlink can be analyzed in the same way. The pro- 
posed model can be justified in the IEEE 802.22 standard, where 
the CR systems are based on cellular basis. By taking advan- 
tage of the parallel OFDMA-based relaying architecture, each 
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Fig. 1 . System architecture of the CRN. 
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Fig. 2. Time slot allocation between a PU and a SU. PU and SU can 
transmit directly (on different subcarriers) or by a cooperation manner. 

PU pair can conduct the bidirectional communication through 
three transmission modes, namely direct transmission, one- and 
two-way relaying, on different sets of subcarriers. As shown 
in Fig. [2 the PUs can transmit directly and the SUs can access 
the PUs' residual subcarriers, or they transmit by a cooperation 
manner. On each cooperated subcarrier, a SU can assist a PU 
(or PU pair) using one- or two-way relaying. This setup can 
fully explore available diversities of the network, including user, 
channel, and transmission mode. 

We model the wireless fading environment by large-scale path 
loss and shadowing, along with small-scale frequency- selective 
fading. The channels between different links experience inde- 
pendent fading and the network operates in slow fading environ- 
ment, so that channel estimation is perfect. We assume that the 
two-hop transmission uses the same subcarrier for both links, 
i.e., the source— >relay link and the relay —^destination link. The 
time slot allocation between a PU and a SU on a cooperated 
subcarrier is illustrated in Fig. [2 in which we further assume 
that the two hops of the cooperative transmission use equal time 
slots. This is true for amplify-and-forward (AF) relaying strat- 
egy because AF needs equal time allocation, but more flexi- 
bility can be provided if the two hops pursue time adaptation 
for decode-and-forward (DF). Nevertheless, we still adopt equal 
time slot allocation between the two hops for simplicity. 

Let N = {1,2,- • -,N} denote the set of subcarriers and 
JC = {!,•■• , K} denote the set of users, with the first 
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Kp being the PU pairs and the remaining Ks = K — Kp 
being SUs. Here k represents PU pair index if 1 < k < 
Kp and represents SU index if Kp + I < k < K. De- 
note k\ and &2 as the two users in the /c-th PU pair, 1 < 
k < K P . Denote P n = [Pi, n , • ■ ',Pk,n, m ' -,^K,n] T and 
R n = [i?i jn , • • •, Rk,m ■ • *j ^K,n] T as the power and achiev- 
able rate vectors on subcarrier n, respectively. If 1 < k < Kp, 
Pk,n = [Pk u n,Pk 2 ,n] T and Rk,n = [Rk u n, Rk 2 ,n] T For in- 
terference avoidance, at most one PU (or PU pair) and one SU 
are active on each subcarrier. This exclusive subcarrier assign- 
ment and best relay selection can be implicitly involved in P n 
and R n . Let P max = [Pf 1 ^, • • P fc max , • • P™ ax ] T denote 
the peak power constraints vector. Again, if 1 < k < Kp, 

pmax = [pmax j pmax^ ^ r = [n, - • ffc, ■ ■ ■, T Kp ] T (with 

^/c = [^fei 7 r fc 2 ] T ) be the rate requirements of the PUs. Denote 
t n = [kpHn,-,^,-,^/ whose elementO < t kjU < 1 
is the duration that SU k transmits on subcarrier n. Note that 
as aforementioned there is at most one SU active on a subcar- 
rier, thus at most one non-zero element in t n . Without loss of 
generality, we assume that additive white noises at all nodes 
are independent circular symmetric complex Gaussian random 
variables, each of them has zero mean and unit variance. As- 
suming channel reciprocity in time-division duplex, we then use 
KmJ 2 > K',bs,„I 2 , and \h p k s uk , n \ 2 (|^ >fc ,,J 2 ) to represent 
the effective channel gains between PU k\ and ki of PU pair k, 
SU k' and BS, and PU k\ fe) and SU k' , respectively, on sub- 
carrier n. For brevity, we denote all of them as a vector H n . A 
PU can cooperate with multiple SUs and a SU can assist multi- 
ple PUs. Thanks to the use of OFDM A, the inter-user interfer- 
ence can be avoided. In addition, the intra-pair interference for 
the PU pairs will be treated as back-propagated self-interference 
and canceled perfectly after two-way relaying. Finally, we let 

P = [Pi, • • "j Pn, ' ' "? Pat] T , R = [Rl, ' ' *, Rn-, ' ' "? Rn] T 

and t = [ti, • • •, t n , • • •, £at] t be the power, achievable rate, and 
time slot allocation matrices, respectively. 

In this paper, our objective is not only to optimally allocate 
power, subcarriers, and time slot but also to choose best trans- 
mission modes and relays for the PUs so as to maximize the 
sum-rate of all SUs while satisfying the individual rate require- 
ment for each of the PUs. Mathematically, the optimization 
problem can be formulated as 
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Rk,n — tk,nC 



^fe,n|^|,BS,n| 2 



,K P + 1 < k < K,Vn, 
(If) 



where C(x) — log 2 (l + x), a| s is the noise variance at the BS, 
and R is the set of achievable rates for the PUs, which is related 
to P n , t n , H n , and the transmission modes. Note again that 
the exclusive subcarrier assignment and best relay selection are 
implicitly involved in P n , R n , and t n . 

Remark 1: In this paper, we assume that a central controller 
is available, so that the network channel state information and 
sensing results can be reliably gathered for centralized process- 
ing. Notice that the centralized CRNs are valid in IEEE 802.22 
standard |[T2ll . where the cognitive systems operate on a cellu- 
lar basis and the central controller can be embedded with a base 
station (BS). This assumption is also reasonable if a spectrum 
broker exists in CRNs for managing spectrum leasing and ac- 
cess E3l[T4) . Such centralized approach is commonly used in a 
variety of CRNs (e.g., ll6 MT3HT9ll ). Compared with distributed 
approaches (e.g., BI4I20I ). a CRN having a central manager that 
possesses detailed information about the wireless network en- 
ables highly efficient network configuration and better enforce- 
ment of a complex set of policies |[T3ll . 

Remark 2: For CRNs, there is no single figure of QoS merit 
to measure the performance of the primary system. In this paper, 
we choose the rate requirement as the QoS metric. Other QoS 
metrics, like outage probability and signal-to-interference-plus- 
noise ratio (SINR), can be easily accommodated in the problem 
formulation. Moreover, the weighted sum-rate maximization for 
the SUs can be taken into account for the fairness issue, which 
does not affect the proposed algorithms in the sequel. 

III. Lagrange Dual Decomposition Based Optimization 

The optimization problem in (Q]) is a mixed integer program- 
ming problem. In ll2T1l . the authors show that for the nonconvex 
resource optimization problems in OFDMA systems, the dual- 
ity gap becomes zero under the time- sharing condition. It is also 
proved in lETll that the time-sharing condition is always satisfied 
as the number of OFDM subcarriers goes to infinity, regardless 
of the nonconvexity of the original problem. This means that 
solving the original problem and solving its dual problem are 
equivalent. Based on the result, the Lagrange dual decomposi- 
tion method is recently applied to OFDMA-based cellular and 
cognitive radio networks in ll22l and fT9lh respectively. In this 
section, we shall apply the result from [CHTl to solve our prob- 
lem in (Q]). In particular, some valuable insights are obtained for 
multi-channel cooperative CRNs, which shows that the general- 
ization from the single-channel case OH) to the multi-channel 
case is nontrivial. 

We first introduce two sets of dual variables, A = [Ai, ■ • 
-,Afc,- • -,\k] t (A*; = [A/ Cl ,A /C2 ] T if 1 < k < K P ) and 
(3= • >, p kr ■ -,Pk p ] t (fa = [/3 kl ^k 2 ] T ) associated with 
constraints ([TBI) and (fTcl) respectively, where A >z and (3 >z 0. 
The Lagrange of the problem in (Q} can be written as 



K N 

L(P,R,t,\,0)= J2 Rk >" 

fe=K P +l n=l 

K / N \ K P / N \ 



k=l 



k=l \n=l 



Define D as the set of all primary variables {P, R, t} that sat- 
isfy constraints (ITdb-dTfb. The dual function is given by 



g(\, (3) = max L(P, R, t, A, /3), 
{p,R,t}eD 

and the dual optimization problem can be expressed as 

min g(\(3) 
s.t. Xh0,(3h0. 

The dual function ([5]) can be rewritten as 



(3) 

(4a) 
(4b) 



TV 



K P 



g(\f3) = $> n (A,/3) + ^ A,P, max - $>* r *' (5) 



/c=l 



/c=l 



where 



g n (\,{3)= max 



_/e=K P + l /e=l 

— ^ ] XkPk,n 
k=l 



(6) 



are the TV independent per-subcarrier-based optimization sub- 
problems. 

Since a dual function is always convex by definition l23ll . 
subgradient-based ellipsoid method ll24l can be used to mini- 
mize #(A, f3) by updating {A, (5} simultaneously along with ap- 
propriate search directions, and it is guaranteed to converge to 
the optimal solution {A*, /3*}. 

Proposition 1: For the dual problem defined in ©, 



TV 



and 



AA fc = P, max 



N 



1 < k < K, 



n=l 



Ap k = Y / R* kin -r k , l<k<K P , 



(7) 



(8) 



are subgradients of g{\ (3), where {P^ n , R^ n } are the optimal 
solutions of © for given {A, /3}. 

Proof: Please see Appendix lAl □ 
As mentioned earlier, there are at most one PU (or PU pair), 
denoted as P, and one SU, denoted as S, active on a subcar- 
rier. Here P and S also represent the best PU (or PU pair) and 
SU respectively, among all possible users, that maximizes © 
for a given subcarrier n. This can be obtained by an exhaus- 
tive search. The complexity is detailed later. Specifically, it 
needs to first compute the optimal powers and rates for all users 
under all transmission modes, then let one PU and/or SU un- 
der one transmission mode that maximizes © active on each 
subcarrier. Therefore, the per- subcarrier problems in © can be 
alternatively expressed as 



max 

P n ,R n ^n 



Rs,n + fipRp,n — X P P Py7l ~ XsPs,< 



S.t. P P , n > 0, P S ,n> 0,0 <t S ,n<l 
Rp,n £ R (Pp,n, Ps,n,ts,n, H n ) 
' BS.nl 2 \ 



Rs,n = ts, n C 



(9a) 

(9b) 
(9c) 

(9d) 



In what follows, for brevity of notation, the subscript n in ^ 
is omitted due to all N per-subcarrier-based subproblems hav- 
ing an identical structure. In addition, for direct transmission 
and one-way relaying, it is observed that the per- subcarrier op- 
timization problem for the two links of a PU pair, i.e., Pi —> P2 
and P2 —> Pi (with or without relaying), has the same structure 
and can be decoupled. Thus, for brevity, we only consider here 
the Pi — > P2 link as an example. In fact, the per- subcarrier op- 
timization for direct transmission and one-way relaying needs 
to first compute the optimal values of the objective function in 
© for both links and then let one of them that has the maximum 
value active on the subcarrier. Moreover, we let 7 = | hJp p 2 n | 2 , 



7i = l^ 1)S ,„l 2 >7 2 = l^ 



S,P 2 ,n\ 



and 7 S 



I^BS^I 2 ' 



A. Direct Transmission 



In this transmission mode, either a PU or a SU occupies solely 
the given subcarrier. The per- subcarrier optimization problem in 
© can be expressed as 

max R s + Pp 2 Rp 2 - X Pl P Pl - X s Ps (10a) 

P Pl >0,P s >0 



s.t. R s = C{P s is) 
Rp 2 =C(P Pi1 ). 



(10b) 
(10c) 



Since the problem in (TTTTb is convex, by applying the Karush- 
Kuhn-Tucker (KKT) conditions l23ll the optimal power alloca- 
tions can be obtained as 



and 



^Pi 



p* 
1 s 



P. 



P2 



aX 



Pi 



\aX s 7s 



(11) 



(12) 



where a = In 2 and = max(0, x). For a given subcarrier, 
the direct transmission further needs to compute the optimal val- 
ues of the objective function in (ITOb over one of Pj^ and Pg with 
the other being zero, and then let one of them that has maximum 
value active. (ITTI) and (Tf2l) show that the optimal power alloca- 
tions are achieved by multi-level water-filling. In particular, the 
water level of each PU depends explicitly on its QoS require- 
ment, and can differ from one another. On the other hand, the 
water levels of all SUs are the same. 

B. One-Way Relaying 

If relaying is required on a given subcarrier, a fraction of 1 — 
ts time is used by a PU to transmit the primary traffic with the 
help of a SU, while the rest ts time is assigned to the SU to 
transmit its own data. In this paper, we focus on DF only for 
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simplicity of presentation, for both one-way relaying and two- 
way relaying. Other relaying strategies are readily applicable to 
our framework and algorithms. The detailed discussion is given 
later. 

In case of DF one-way relaying, the per-subcarrier problem 
in © can be rewritten as 



max R s + (3p 2 Rp 2 - ^p 1 Pp 1 - >^sPs 

P Pl ,Ps,ts 

s.t. P Pl > 0, P s > 0, < t s < 1 
Rs = tsC(Psjs) 
(l-*s) 



(13a) 

(13b) 
(13c) 



Rp 2 



min {C(P Pl7 i), C(P Pll + P572)} 

(13d) 



In (I13dlh the first term in the min-operation is the achievable 
rate of the P\ —> S link, and the second term is the achievable 
rate by maximum ratio combining between the S — >• P2 link 
and P\ —> P2 link. The following proposition is established for 
determining the optimal value of the time slot allocation variable 
ts. 

Proposition 2: For each cooperated subcarrier, a SU exclu- 
sively acts as a relay for cooperative transmission or transmits 
traffic for itself. 

Proof: The proposition means that the time slot allocation 
variable ts is binary, i.e., t$ G {0, 1}, which can be proved by 
contradiction. 

Assume that the optimal solution of (13]) is (t* s , P Pi , Pg) with 
< tg < 1. Next, we show that we can always find another 
better solution with ts being binary. 

Let us rewrite the objective function d!3ab as 



f(t* s ,P* Pl ,P£) = t s 



-p P2 min |^C(P^ 7 i), \c{P* Pll + P* sl2 ) 
+p P2 min |ic(P^ l7 i), \c{P* Pll + P* sl 2) 

-\ Pl P Pl - X S Pg. 



(14) 



lfC(P* s7s ) < pp 2 mm{\C{P Pill ),\C{P* Pil + P* sl2 )}, we 
have 



\ Pl P Pl -\sP*s 

Pp 2 



This proposition also holds for two-way relaying as discussed 
in the next subsection. The proof is similar and hence ignored. 

Proposition [2] significantly simplifies the per-subcarrier opti- 
mization problem in (IT3l) without loss of optimality by an ex- 
haustive search over ts. Specifically, we set ts = and ts = 1 
to compute the optimal values of (113aK respectively, then follow 
the one that has the maximum value. 

The intuition is that, on a cooperated subcarrier (see Fig. [2]), 
if the subcarrier condition on the SU— ^BS link is good but on 
the cooperative link is poor, it is better that the PU leases the 
whole transmission time slot to the SU. Otherwise, the SU com- 
pletely devotes itself as a relay to the PU. In other words, if a 
SU exclusively forwards a PU's traffic on a subcarrier, the PU 
shall lease other subcarrier(s) to the SU as remuneration. This 
channel-swap based multi-channel cooperation fundamentally 
differs from the single-channel cooperation case 0141 where if 
a SU forwards a PU's traffic, it must benefit from the PU on 
the channel. The spectral efficiency improvement brings more 
cooperation opportunities and leased subcarriers, and thus, the 
total throughput of the secondary system is increased. In the 
following we consider ts = and ts = 1, respectively. 

Case 1: ts = 0. In this case, S exclusively acts as a relay on 
a cooperated subcarrier. In DF one-way relaying, it is intuitive 
that Rp 2 is maximized when C (P Pl 71) = C(Pp 1 j + ^572), 
which leads to 

Ps = iP Pl , (17) 

where 7' = (71 — 7V72. It is noted that DF occurs only if 
71 > 7. Substituting ([171) to (IT3l) and let ts = 0, the problem 
can be rewritten as 

max (3p 2 Rp 2 - (A Pl + Xsl)Pp 1 (18a) 
Pp 1 >o 



s.t. R P2 = -C(P Plll ). 



(18b) 



The above is a convex problem. By applying the KKT condi- 
tions, the optimal power allocation is given by 



Pi 



Pp 2 



1 



2a(X Pl + A57') 71 



(19) 



(15) 



Similarly, if C(P^ S ) > ^ min {C(P* i7 i), C(P^ 7 + ^72)} 
we have 

f(thP Pl :Ps) < C(P* sl * s )-\ Pl P^-\sP*s 
< C(P* sl * s ) - \ s P*s 
= /(1,0,P5). (16) 

These results contradict the assumption. This completes the 
proof. □ 



Pg can be obtained according to ([T71) . The above optimal power 
allocation (|T9t shows that higher channel gain 71, meaning 
higher 7', results in lower water level, which is the extra fea- 
ture compared with the standard water-filling approach (e.g., 
(TTTT) and (fT2l) in direct transmission). One also observes that 
lower channel gain 72 leads to lower water level and vice versa. 

Case 2: ts = 1. In this case, S uses a cooperated subcarrier 
solely for its own transmission. The optimal power allocation 
can be easily obtained and is the same as (fT2l) . 

C Two-Way Relaying 

The two-way communication between Pi and P2 assisted by 
S takes place in two phases. Specifically, in the first phase, also 
known as multiple-access (MAC) phase, Pi and P2 concurrently 
transmit signals to the assisting S. In the second phase, known 
as broadcast (BC) phase, S broadcasts the processed signals to 
both Pi and P2 . Different from direct transmission and one-way 
relaying, two-way relaying must occur in pair ll25l - [30l . Thus 
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both the bidirectional links are taken into account together for 
two-way relaying. Here we only analyze the case of ts = 0, and 
the case of ts = 1 is omitted since the optimal power allocation 
is the same as (IT2l) if ts = 1. 

The per-subcarrier problem in © can be expressed as (recall 
that t s 0) 

max Bp, Rp, + 0p 9 Rp 9 

Pp>0,P s >0,R P 

- X Pl P Pl - Xp 2 Pp 2 - X S P S (20a) 

s.t. R P G R (Pp , Ps , 71 ? 72) = CmacH^bc, 

(20b) 

where Cmac and Cbc are the capacity regions for the MAC and 
BC phases, respectively (T2EUE2I • Specifically, 

Cmac = {[R Pl Rp 2 ]\r Pi < \c(Pp 2l2 ),Rp 2 < \c(P Plll ), 
R Pi +Rp 2 < \c(Pp 2 72 + Pp, 71)}, (21) 

and 



BC 



{[P Pl Pp 2 ]|p Pi < \c(P sll ),R P2 < l -C(P sl2 )). 

(22) 

Note that the channel reciprocity is used in the BC phase, which 
is justified by the time-division duplex mode. Since both Cmac 
and Cbc are convex sets, and so is their intersection, the prob- 
lem in (l20l) is a convex problem and can be solved by convex 
techniques. 

Let ql\ and a 2 be the two dual variables associated with the 
two rate constraints in d22l) . We first incorporate the two rate 
constraints in (l22b into the objective function and rewrite the 
Lagrange dual problem of (f20l) as 

min max /3p 1 Rp 1 + /3p 2 Rp 2 

on,a 2 {Pp,Rp}eC DF 

-X Pi P Pi -Xp 2 Pp 2 -X s P s 



Rpl - \c(Psli) 

1 



Rp, - ^c(p sl2 ) 



s.t. 



— a 2 
ai >0,a 2 > 0, 



(23a) 
(23b) 



where Cdf is the set of the remaining constraints in (l20l) that 
{Pp, Pp} must satisfy. The minimization over {a± : a 2 } can be 
done using ellipsoid method with the fact that ^C(Psji) — Rp 1 
and \C(Ps~f 2 ) — Rp 2 are subgradients of a\ and a 2 , respec- 
tively. It is observed that the optimization variables in (l23l) are 
separable. Therefore, the maximization over {Pp, Pp} in (f23l) 
can be decomposed into two subproblems that can be solved 
separately. The two subproblems are 

max (/3 Pl - a 1 )R Pl + (/?p 2 - a 2 )^p 2 

-A Pl P Pl -A P2 Pp 2 (24a) 
s.t. {Pp,i? P } G Cmac, (24b) 



and 

max ^C(Ps7i) + ^C(Ps7 2 )-\sPs- (25) 

For brevity of notation, we let a[ = (3 Pl — ql\ and a' 2 = 
f3p 2 — a 2 . It is noted that both a[ and a' 2 must be nonnega- 
tive, i.e., f3p 1 > ol\ and f3p 2 > a 2 . In the following we present 
the solution to each subproblem. 

The subproblem in (l24l) is a classic resource allocation prob- 
lem in the Gaussian MAC (33), where the optimal power and 
rate allocations can be achieved by successive decoding. Specif- 
ically, users' signals are decoded one by one in an increasing rate 
weight order l33l . Without loss of generality, we assume that 
ot'i > ol 2 (here a[ and a' 2 can be regarded as the rate weights 
for Pi and P 2 , respectively). Based on the polymatroid struc- 
ture of the Gaussian MAC ll33lh we then incorporate the three 
rate constraints in (fzTI) into the objective function of (l24lh the 
subproblem in (f24l) can be expressed as 

max y^(Pp 2 72) + ^ [C(P P2 72 + Pp i7 i) - C(Pp 2 72)] 



-Ap.Pp, 



Ap 2 Pp 2 



(26) 



It is easy to validate that the objective function of (l26l) is jointly 
concave in Pp 1 and Pp 2 . By applying the KKT conditions, the 
optimal power allocations can be obtained as 



and 



p* 

^P<2 



^Pl - 



K - 0*2)71 _ J_ 

_2a(7iAp 2 - 72APJ 72 



1 + 



1 

2a 



a f 2 



(27) 



(28) 



_Xp 1 7iAp 2 — 72Ap x _ 

It is observed that the optimal power allocations in the MAC 
phase have the form of water- filling. Moreover, the follow- 
ing proposition is provided according to the above closed-form 
power allocations. 

Proposition 3: For a[ > a 2 , a necessary condition for the 
occurrence of two-way relaying is 71 Ap 2 > 72 X Pl . 

Proof: Please see Appendix |B] □ 

The subproblem in (l25t is also convex since its objective 
function is concave in Ps. By applying the KKT conditions, 
the optimal power allocation is given by 



-46>i 



20i 



0, 



if As < 
otherwise 



2a , (29) 



where 6\ = 2aXsjij 2 , 2 = 2aXs(ji + 72) - 7172(^1 + ^2), 
and 63 = 2aXs — o^iji — a 2 j 2 . 

Remark 3: For both one- and two-way relaying, direct trans- 
mission mode is optimal if P| =0. In this case, the optimal 
power and rate allocations are the same as those obtained in the 
direct transmission mode in Section III-A. Moreover, it is noted 
that two-way relaying occurs if P| , Pp , and Pp 2 are all posi- 
tive. For two-way relaying, another interesting case is that P^ is 
positive and one of Pp and Pp 2 is equal to zero. In other words, 
one direction is inactive. In this case, one-way relaying must be 
optimal. 

Remark 4: In our centralized framework, the cooperation is 
between the primary system and secondary system rather than 
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among individual users. Thus some SUs may not transmit their 
own traffic because they are the best option for primary traffic 
relaying. In this case, these SUs may be re-scheduled for trans- 
mission at the next transmission frame by a higher layer sched- 
uler for long-term fairness. However, analysis on higher layer 
scheduling is beyond of the scope of this paper. 

Remark 5: In this paper, we employ DF just for an illustra- 
tion purpose. Other relaying strategies, like AF and compress- 
and-forward (CF), are generally applicable to our proposed 
framework and algorithms. However, the achievable rate ex- 
pressions of AF and CF are not concave, for both one- and two- 
way relaying @ To overcome the difficulty, some approximations 
can be adopted for AF and CF such that the achievable rate ex- 
pressions are concave, and thus they can be solved using convex 
optimization techniques. 

After obtaining the optimal solution in the dual domain, we 
now need to obtain the optimal solution to the original primal 
problem in (Q]). Due to the non-zero duality gap, the optimal 
solution obtained in the dual domain may not satisfy all the 
constraints in the original primal problem. To tackle this prob- 
lem, we first obtain the optimal transmission mode selection and 
user-assignment for each subcarrier using the method in the dual 
domain, then the primal problem in (Q]) reduces to a pure power 
allocation problem and it is convex. By applying KKT con- 
ditions, the optimal power allocations follow the same expres- 
sions in the dual domain and the details are omitted here. This 
approach is asymptotically optimal due to the vanishing duality 
gap when the number of subcarriers is sufficiently large ll2T1l . 

At the end of the section, we analyze the computational 
complexity of the proposed algorithm. The complexity of de- 
termining P and S on each subcarrier for direct transmis- 
sion is 0(2Kp + Ks), and for one- and two-way relaying 
are 0(2K P K S ) and 0(K P K s ), respectively. Note that the 
complexity of the search over ts = 1 is implicitly contained 
in the optimization of direct transmission mode. Therefore, 
the total complexity of solving all TV per- subcarrier problems 
is O (N(2K P + K s + 3K P K S )). Combining the complex- 
ity of the ellipsoid method, the total complexity of solving the 
dual problem is O (N(2K P + K s + 3K P K S )(4K P + K s ) 2 ), 
which is linear in N and polynomial in K P and Ks. 

IV. Simulation Results 

In this section, we evaluate the performance of the pro- 
posed cooperative scheme using simulation. The conven- 
tional scheme without cooperation is selected as a bench- 
mark, which corresponds to the optimization of direct trans- 
mission in our proposed algorithm and its complexity is 
O (N(2K P + K S )(4K P + K s ) 2 ). As another benchmark, the 
performance of the Fixed Transmission Mode (FTM) based al- 
location is also presented. In particular, the FTM scheme lets 
the transmission mode for each PU be pre-fixed according to 
nodes' geographical information, and other optimizations are 
the same with the proposed optimal algorithm. This is attractive 
for practical systems where path loss dominates the performance 
of the network nodes. In specific, a PU is assigned to the direct 
transmission mode if the path loss (or distance) of the source- 

2 It is noted that the nonconvexity does not affect Proposition [2] 
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Fig. 3. Sum-rate of the secondary system versus transmit SNR per 
subcarrier, with K P = 2 PU pairs, K s = 4 SUs and rate requirement 5 
bit/OFDM symbol for all PUs. 

destination link is smaller than that of all source-relay links and 
the cooperative transmission modes are used otherwise. When 
it is assigned the cooperative transmission modes, two-way re- 
laying is adopted if the path losses of the source-relay link and 
the relay-destination link is about the same, otherwise one-way 
relaying is used (in this case, Pi — >> P2 direction is performed). 
Note that if a PU is assigned two-way relaying, the other PU in 
the same pair is also assigned two-way relaying. For those PUs 
who need SUs' assistance, we assign the nearest SU to each 
PU, the search over the suitable SU for each PU is reduced to 
0(1). Hence, the total complexity of this suboptimal algorithm 
is O (N(5K P + K S )(4K P + K s ) 2 ). 

We consider an a primary network in a 1 km by 1 km square 
area and a cellular secondary network whose BS is located in 
the center of the square and with 1 km radius. All users are 
randomly but uniformly distributed. The statistical path loss 
model and shadowing are referred to l34l . where we set the path 
loss exponent to be 4 and the standard deviation of log-normal 
shadowing to be 5.8 dB. The small-scale fading is modeled by 
Rayleigh fading process, where the power delay profile is ex- 
ponentially decaying with maximum delay spread of 5 /is. A 
total of 2000 independent channel realizations were used. Dif- 
ferent channel realizations are with different node locations. We 
set the number of OFDM subcarriers be N = 64. Without loss 
of generality, we let all users have the same maximum power 
constraints, and all PUs have the same rate requirements. In all 
of the simulations, we fix K P = 2 PU pairs (or equivalently 4 
PUs) in the network. Without loss of generality, we let all users 
have the same peak power constraints. 

Figs. [3] and HI compare the sum-rate of the secondary system 
versus user peak power (in dB) achieved by different schemes 
when there are Ks = 4 and Ks = 8 SUs, respectively. In both 
figures, the PU rate requirement is 5 bits/OFDM symbol for each 
PU, and the dual optimum values serve as the performance up- 
per bounds. It is first observed that the proposed cooperation 
scheme approaches the upper bound very tightly, which veri- 
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Fig. 4. Sum-rate of the secondary system versus peak power per 
subcarrier, with K P = 2 PU pairs, K s = 8 SUs and rate requirement 5 
bit/OFDM symbol for all PUs. 

fies the effectiveness of the proposed algorithm. One also ob- 
serves that the proposed scheme outperforms the conventional 
non-cooperative scheme by a significant margin. In particular, 
compared with the conventional scheme, about 60% throughput 
improvement is achieved in our proposed scheme. The tremen- 
dous improvement is as the remuneration for cooperative diver- 
sity, selection diversity, and network coding gain that the sec- 
ondary systems provides to the primary system. Second, one 
also observes that our proposed scheme improves 20% through- 
put over the FTM scheme. This clearly suggests the benefits of 
bidirectional transmission mode adaptation and SU selection for 
the PUs. Third, from Fig. [4] with Ks = 8 SUs, it is observed 
that our proposed scheme also outperforms the FTM and con- 
ventional schemes substantially. Note that a larger Ks results in 
higher computational complexity mainly due to the updates of 
dual variables. 

We next study the sum-rate of the secondary system versus 
the different PU rate requirements in Fig. where we fix trans- 
mit SNR 10 dB and Ks = 4 SUs. As expected, our proposed 
scheme outperforms the FTM and conventional schemes con- 
siderably over all ranges of PU rate requirements. This further 
demonstrates the effectiveness of the proposed scheme. 

V. Conclusion 

This paper studied the OFDMA-based bidirectional CRNs 
with cooperation between the primary and secondary systems, 
for supporting communication services with diverse QoS re- 
quirements. We proposed an optimization framework for joint 
optimization of bidirectional transmission mode selection, SU 
selection, subcarrier assignment, power control, and time slot al- 
location. We converted this mix integer programming problem 
with exponential complexity into a convex problem using the 
dual decomposition method and developed efficient algorithms 
with polynomial complexity. 

A few important conclusions have been made throughout this 
paper. Firstly, the time slot allocation between a PU and a 



Fig. 5. Sum-rate of the secondary system versus PU rate requirement, 
with K P = 2 PU pairs, K s = 4 SUs and peak power is 10 dB per 
subcarrier. 

SU on a cooperated subcarrier is binary. Secondly, the pro- 
posed framework can greatly improve the total throughput of 
the secondary system by about 60%, compared with the non- 
cooperative scheme. Thirdly, choosing the appropriate trans- 
mission modes for the PUs is necessary. Last but not least, 
transmission mode adaptation and SU selection over different 
subcarriers can enhance the total throughput by about 20%. 

The proposed algorithm can be used as the performance upper 
bound for suboptimal or distributed algorithms. In future work, 
it will be interesting to investigate incentive-based distributed 
schemes. 

APPENDICES 
I. Proof of Proposition Q] 

By definition of g( A, (3) in ([3]), we have 

K N K / N \ 

5 (v, 0)> Yl E R in + E a m p " ax - E P ln 

fc=K P +ln=l k=l \ n=l J 

Kp I N \ 

+ \J2 R k,n-rk) 

fc=l \n=l ) 

K / N \ 

= g(X,0) + £(A£ - A fc ) Pf ax - £ P fc *„ 

k=l \ n=l J 

K P / N \ 

+ I>*-fl0 [Y, R in-r k . (30) 



k=l 



Hence, Proposition \T\ is proven by using the definition of sub- 
gradient. 

II. Proof of Proposition [3] 

Both Pp 2 and Pp 2 must be positive if two-way relaying oc- 
curs, besides Pg is positive. We first investigate Pp 2 in (l27l) . It 
is easy to observe that 7iAp 2 must be greater than 72Ap 1? oth- 
erwise the first term in (l27t is negative, and thus, Pp = 0. For 
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Pp 2 in (l28k we let the first term is greater than that of the sec- 
ond term, i.e., a' 2 /\p 1 > (a[ — 0^2)72/(71 Ap 2 — 72ApJ. After 
some manipulations, we obtain / yi\p 2 /j2^p 1 > a[/a 2 - Com- 
bining the condition a[/a f 2 > 1, we obtain 7iAp 2 > 72 Ap^ 
This completes the proof. 
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